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Abstract. An accurate analytic approximation of the 
transfer function for the power spectra of primordial den- 
sity perturbations in mixed dark matter models is pre- 
sented. The fitting formula in a matter-dominated Uni- 
verse (Sl = Qm = 1) is a function of wavenumber k, 
redshift z and four cosmological parameters: the density 
of massive neutrinos, Q u , the number of massive neutrino 
species, N v , the baryon density, and the dimensionless 
Hubble constant, h. Our formula is accurate in a broad 
range of parameters: k < 100 h/Mpc, z < 30, Q, v < 0.5, 
N v < 3, tt b < 0.3, 0.3 < h < 0.7. The deviation of the 
variance of density fluctuations calculated with our for- 
mula from numerical results obtained with CMBfast is less 
than 6% for the entire range of parameters. It increases 
with i} b h 2 and is less than < 3% for Q b h 2 < 0.05. 

The performance of the analytic approximation 
of MDM power spectra proposed here is com- 
pared with other approximations found in the lit- 
erature (Holtzman 1989, Pogosyan & Starobinsky 1995, 
Ma 1996, Eisenstein & Hu 1997b). Our approximation 
turns out to be closest to numerical results in the pa- 
rameter space considered here. 
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1. Introduction 

Finding a viable model for the formation of large scale 
structure (LSS) is an important problem in cosmology. 
Models with a minimal number of free parameters, such 
as standard cold dark matter (sCDM) or standard cold 
plus hot, mixed dark matter (sMDM) only marginally 
match observational data. Better agreement between pre- 
dictions and observational data can be achieved in models 
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with a larger numbers of parameters (CDM or MDM with 
baryons, tilt of primordial power spectrum, 3D curvature, 
cosmological constant, see, e.g., Valdarnini et al. 1998 and 
refs. therein). In view of the growing amount of obser- 
vational data, we seriously have to discuss the precise 
quantitative differences between theory and observations 
for the whole class of available models by varying all 
the input parameters such as the tilt of primordial spec- 
trum, n, the density of cold dark matter, £Icdm, hot 
dark matter, f2„, and baryons, f2b, the vacuum energy or 
cosmological constant, and the Hubble parameter h 
(h = i/o/100 km/s/Mpc), to find the values which agree 
best with observations of large scale structure (or even to 
exclude the whole family of models.). 

Publicly available fast codes to calculate the 
transfer function and power spectrum of fluctua- 
tions in the cosmic microwave background (CMB) 
(Seljak & Zaldarriaga 1996, CMBfast) are an essential in- 
gredient in this process. But even CMBfast is too bulky 
and too slow for an effective search of cosmological pa- 
rameters by means of a ^-minimization, like that of Mar- 
quardt (see Press et al. 1992 ). To solve this problem, an- 
alytic approximations of the transfer function are of great 
value. Recently, such an approximation has been pro- 
posed by Eisenstein & Hu 1997b (this reference is denoted 
by EH2 in the sequel). Previously, approximations by 
Holtzman 1989, Pogosyan & Starobinsky 1995, Ma 1996 
have been used. 

Holtzman's approximation is very accurate but it 
is an approximation for fixed cosmological parame- 
ters. Therefore it can not be useful for the pur- 
pose mentioned above. The analytic approximation 
by Pogosyan & Starobinsky 1995 is valid in the 2- 
dimensional parameter space h), and z (the redshift). 
It has the correct asymptotic behavior at small and large 
k, but the systematic error of the transfer function T(k) is 
relatively large (10%-15%) in the important range of scales 
0.3 < k < 10 /i/Mpc. This error, however introduces dis- 
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crepancies of 4% to 10% in ctr which represents an inte- 
gral over k. Ma's analytic approximation is slightly more 
accurate in this range, but has an incorrect asymptotic 
behavior at large k, hence it cannot be used for the analy- 
sis of the formation of small scale objects (QSO, damped 
Ly a systems, Ly a clouds etc.). 

Another weak point of these analytic approxima- 
tions is their lack of dependence on the baryon den- 
sity. Sugiyama's correction of the CDM transfer func- 
tion in the presence of baryons (Bardeen et al. 1986, 
Sugiyama 1995) works well only for low baryonic con- 
tent. Recent data on the high-redshift deuterium abun- 
dance (Tytler ct al. 1996), on clustering at lOOMpc/ft. 
(Eisenstein et al. 1997) and new theoretical interpreta- 
tions of the Ly a forest (Weinberg et al. 1997) suggest that 
Sib may be higher than the standard nucleosynthesis value. 
Therefore pure CDM and MDM models have to be mod- 
ified. (Instead of raising fib, one can also look for other 
solutions, like, e.g. a cosmological constant, see below.) 

For CDM this has been achieved by Eisenstein & Hu 
(1996, 1997a 1 ) using an analytical approach for the de- 
scription of small scale cosmological perturbations in the 
photon-baryon-CDM system. Their analytic approxima- 
tion for the matter transfer function in 2-dimensional pa- 
rameter space (r^M^ 2 , ^b/^Af) reproduces acoustic oscil- 
lations, and is quite accurate for z < 30 (the residuals 
are smaller than 5%) in the range 0.025 < f^M^ 2 < 0.25, 

< fib/ & m < 0.5, where Qm is the matter density pa- 
rameter. 

In EH2 an analytic approximation of the mat- 
ter transfer function for MDM models is proposed 
for a wide range of parameters (0.06 < 51m^ 2 < 

o.4, n b /n M < o.3, n u /n M < 0.3 and z < 30). 

It is more accurate than previous approximations by 
Pogosyan & Starobinsky 1995, Ma 1996 but not as pre- 
cise as the one for the CDM+baryon model. The baryon 
oscillations are mimicked by a smooth function, therefore 
the approximation looses accuracy in the important range 
0.03 < k < 0.5 h/Mpc. For the parameter choice VLm = 1, 
VL U = 0.2, SI;, = 0.12, h = 0.5, e.g., the systematic resid- 
uals are about 6% on these scales. For higher f2„ and f2& 
they become even larger. 

For models with cosmological constant, the motiva- 
tion to go to high values for Q u and Qb is lost, and the 
parameter space investigated in EH2 is sufficient. Models 
without cosmological constant, however, tend to require 
relatively high baryon or HDM content. In this paper, 
our goal is thus to construct a very precise analytic ap- 
proximation for the redshift dependent transfer function 
in the 4-dimcnsional space of spatially flat matter domi- 
nated MDM models, Tmdm^; f2„, N u , fib, h; z), which is 
valid for SIm = 1 and allows for high values of f2„ and fit,. 
In order to keep the baryonic features, we will use the EH1 
transfer function for the cold particles+baryon system, 

1 This reference is denoted by EH1 in this paper. 



TcDM+b(k; Qb,h), and then correct it for the presence of 
HDM by a function D(k\ Q„, N„, O;,, h; z), making use of 
the exact asymptotic solutions. The resulting MDM trans- 
fer function is the product T MDM (k) — T C DM+b(k)D(k). 

To compare our approximation with the numerical re- 
sult, we use the publicly available code 'CMBfast' by Sel- 
jak & Zaldarriaga 1996. 

The paper is organized as follows: In Section 2 a short 
description of the physical parameters which affect the 
shape of the MDM transfer function is given. In Section 
3 we derive the analytic approximation for the function 
D(k). The precision of our approximation for TMDAi(k), 
the parameter range where it is applicable, and a compar- 
ison with the other results are discussed in Sections 4 and 
5. In Section 6 we present our conclusions. 

2. Physical scales which determine the form of 
MDM transfer function 

We assume the usual cosmological paradigm: scalar pri- 
mordial density perturbations which are generated in the 
early Universe, evolve in a multicomponent medium of 
relativistic (photons and massless neutrinos) and non- 
relativistic (baryons, massive neutrinos and CDM) parti- 
cles. Non-relativistic matter dominates the density today, 
Q.M = ftb + ftf + Sic dm- This model is usually called 
'mixed dark matter' (MDM). The total energy density 
may also include a vacuum energy, so that U = SIm + ^A- 
However, for reasons mentioned in the introduction, here 
we investigate the case of a matter-dominated flat Uni- 
verse with SIm = 1 and = 0. Even though Oa 
seems to be favored by some of the present data, our main 
point, allowing for high values of fij,, is not important in 
this case and the approximations by EH2 can be used. 

Models with hot dark matter or MDM have been 
described in the literature by Fang, Xiang & Li 1984, 
Shafi & Stecker 1984, Valdarnini & Bonometto 1985, 
Holtzman 1989, Lukash 1991, Davis, Summers & Schlegel 
1992, Schaefer & Shafi 1992, Van Dalen & Schaefer 1992, 
Pogosyan & Starobinsky 1993, 1995, Novosyadlyj 1994, 
Ma & Bertschinger 1994, 1995, Seljak & Zaldarriaga 1996, 
EH2, Valdarnini et al. 1998 and refs. therein. Below, we 
simply present the physical parameters which determine 
the shape of the MDM transfer function and which will 
be used explicitly in the approximation which we derive 
here 2 . 

2 Recall the definitions and relationship between the MDM 
and the partial transfer functions 

Tmdm = QcdmTcdm + S~lvT u + fi^Ti, , 

1(K) - 8(0,z)5(k,z ln ) ' 

where S(k, z) is the density perturbations in a given component 
and z in is a very high redshift at which all scales of interest 
are still super horizon. 
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Fig. 1. The transfer function of density perturbations of CDM, 
baryons and HDM at z = 10 (calculated numerically). 



Since cosmological perturbations cannot grow signif- 
icantly in a radiation dominated universe, an important 
parameter is the time of equality between the densities of 
matter and radiation 



2.4 x 10 4 2 4 
1-^/7.4 ^ 
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(1) 



where i 7 = T 1 /2.726K is the CMB temperature today, 
N v =l, 2 or 3 is the number of species of massive neutrinos 
with equal mass (the number of massless neutrino species 
is then 3 — N v ). The scale of the particle horizon at this 
epoch, 

(2) 



keq 



4.7 x 10 _4 v /l + z e9 h/Mpc, 



is imprinted in the matter transfer function: perturbations 
on smaller scales (k > k eq ) can only start growing after 
z eq , while those on larger scales (k < k eq ) keep growing 
at any time. This leads to the suppression of the transfer 
function at k > k eq . After z eq the fluctuations in the CDM 
component are gravitationally unstable on all scales. The 
scale k eq is thus the single physical parameter which de- 
termines the form of the CDM transfer function. 

The transfer function for HDM (y) is more compli- 
cated because two more physical scales enter the problem. 
The time and horizon scale when neutrino become non- 
relativistic [m v ~ 3T„) are given by 



Fig. 2. The same as in Fig.l but for z=0. 



k n 



3.3 x 10" 



y^l + avXl + Zeg) h/Mpc, (3) 

where x v = Ct v /Q veq , VL W eq ~ N„/(7A — N„) is 
the density parameter for a neutrino component becom- 
ing non-relativistic just at z eq . The neutrino mass can 
be expressed in terms of Sl v and N u as (Peebles 1993) 
m v = 94fl v h 2 N-H~ 3 eV. 

The neutrino free-streaming (or Jeans 3 ) scale at z < 



is 



k F (z) ~ 59. 



l + z e 



l + z 



Q„N-H- 4 h 3 /Mpc, (4) 



which corresponds to the distance a neutrino travels in 
one Hubble time, with the characteristic velocity v v ~ 



ar i and k nr ^k eq for fl^^Oj, 



eg 



The amplitude of ^-density perturbation on small 
scales (k > k nr ) is reduced in comparison with large scales 
(k < k nr ). For scales larger than the free-streaming scale 
(k < hp) the amplitude of density perturbations grows in 
all components like (1 + z)' 1 after z eq . Perturbations on 
scales below the free-streaming scale (k > k F ) are sup- 
pressed by free streaming which is imprinted in the trans- 
fer function of HDM. Thus the latter should be parame- 
terized by two ratios: k/k nr and k/kp. 

3 Formally the Jeans scale is 22.5% less than the free- 
streaming scale (Bond & Szalay 1983, Davis, Summers & 
Schlegel 1992), however, kF is the relevant physical parame- 
ter for collisionless neutrini. 
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The transfer function of the baryon component is de- 
termined by the sound horizon and the Silk damping scale 
at the time of recombination (for details see EH1). 

In reality the transfer function of each component is 
more complicated due to interactions between them. At 
late time (z < 20), the baryonic transfer function is 
practically equal to the one of CDM, for models with 
fib < flcDM (see Figs. 1,2). After z eq , the free-streaming 
scale decreases with time (neutrino momenta decay with 
the expansion of the Universe whereas the Hubble time 
grows only as the square root of the scale factor, see 
Eq. (4)), and neutrino density perturbations at smaller 
and smaller scales become gravitationally unstable and 
cluster with the CDM+baryon component. Today the v 
free-streaming scale may lie in the range of galaxy to clus- 
ters scales depending on the v mass. On smaller scales the 
growing mode of perturbation is concentrated in the CDM 
and baryon components. Matter density perturbation on 
these scales grow like ~ t a , where a = (\/25 — 24f2„ — 1)/6 
(Doroshkevich et al. 1980). 

3. An analytic approximation for the MDM trans- 
fer function 

To construct the MDM transfer function we use the an- 
alytic approximation of EH1 for the transfer function of 
cold particles plus baryons and correct it for the presence 
of a ^-component like Pogosyan & Starobinsky 1995 and 
Ma 1996: 

Tmdm (k) = T CD M+b{k)D{k) . (5) 
The function D(k) must have the following asymptotics: 
D(k < k nr ) = 1, 





D{k^>k F ) = {l-VL v ) 



l + z 

1 



where (3 = 1 — 1.5a. After some numerical experimenta- 
tion we arrive at the following ansatz which satisfies these 
asymptotics 



D(k) = 



i + (i - a,) 1 //^ 



1 + {a 4 k/k nr ) 3 



(6) 

We minimize the residuals in intermediate region 
(k nr < k < kp) by determining ca as best fit coefficients 
by comparison with the numerical results. 

By x 2 minimization (Press et al. 1992 ) we first deter- 
mine the dependence of the coefficients on on VL V keeping 
all other parameters fixed, to obtain an analytic approx- 
imation a.i(Q. v ,z). The main dependence of Tmdm {k) on 
fii,, N v , h and z is taken care of by the dependence of 
TcDM+b, k nr , &f and of the asymptotic solution on these 
parameters. We then correct a, by minimization of the 



Q 



0.1 



1 1 rTTTTTj 1 1 — 1 1 1 1 1 1 1 1 1 


1 1 1 1 ll| 1 - 


-i — i i 1 1 ii| 1 — i — n 






Ov=0.1 . 


- 




\. 0.2 ■ 






\. 0.3 


_ 




\ 0.4 ' 


o.b=0.06, h=0.5, z=0 






\ 0.5 




1E-3 



0.01 



0.1 1 

k (h/Mpc) 



10 



Fig. 3. D(k) = TMDM(k)/TcDM+b(k) as calculated numeri- 
cally (solid line) and our analytic approximation (dotted line) 
for different values of Q. v . The numerical results and the ap- 
proximations overlay perfectly. 



residuals to include the slight dependence on these pa- 
rameters. 

Finally, the correction coefficients have the following 
form: 

a, = aiA i {z)B i {n h )Ci{h)D i {N v ), 

where a, = a 4 (f>„), ,4,(0) = B 4 (0.06) = C 4 (0.5) - D t (l) = 
1. The functions Ai depend also on Q„. 

For all our calculations we assume a CMB temperature 
of T 7 = 2.726if (Mather et al. 1994, Kogut et al. 1996). 

3.1. Dependence on fl u and z. 

We first set h = 0.5, f2& = 0.06, N u = 1 and determine 
for Q. v = 0.1, 0.2, 0.3, 0.4, 0.5 and z = 0, 10, 20. We 
then approximate D(k) by setting on = CLiAi(z), where 
Ai(z) = (1 + z ) b > +c * (1+z l The dependences of a t , b t and 
Cj on Q„, as well as Bi(Qb), Ci(h) and Di{N v ) are given 
in the Appendix. The functions D(k) for different f2„ and 
its fractional residuals are shown in Figs. 3 and 4. 

We now analyze the accuracy of our analytic ap- 
proximation for the MDM transfer function Tmdm(^) = 
TcDM+bD(k) which in addition to the errors in D{k) 
contains also those of TcDM+b(k) (EH1). We de- 
fine the fractional residuals for Tmdm(^) by (T(k) — 
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Fig. 4. Fractional residuals of D(k) given Fig. 5. The numerical MDM transfer function 
by (-D(fc) — DcMBfast(k))/DcMBfast(k) for different values at z = (thick solid lines) and the analytic approximations 



of f2„. 



TMDAi(k) = TcDM+bD(k) (dotted thin lines, which perfectly 
overlay with the numerical result). 



TcMBfast{k))/T C MBfast(k). In Fig. 5 the numerical re- 
sult for TMDHi(k) (thick solid lines) and the analytic ap- 
proximations (dotted thin lines) are shown for different 
f2„. The fractional residuals for the latter are given in 
Fig. 6. Our analytic approximation of TmdmCO is suffi- 
ciently accurate for a wide range of redshifts (see Fig. 7). 
For z < 30 the fractional residuals do not change by more 
than 2% and stay small. 

3.2. Dependence on fi& and h. 

We now vary fib fixing different values of Q, v and setting 
the other parameters h = 0.5, N v = 1. We analyze the 
ratio D(k;Cl v ,Clb)/D(k;Cl v ,Clb = 0.06). Since the domi- 
nant dependence of Tmdm(^) on fij, is already taken care 
of in TcDM+b(k), D(k) is only slightly corrected for this 
parameter. Correction factors Bi(Qb) (~ 1) as a second 
order polynomial dependence on ilt,/0.06 with best-fit co- 
efficients are presented in the Appendix. The fractional 
residuals of TMDAf(k) for different f2& are shown in Fig. 8. 

The maximum of the residuals grows for higher baryon 
fractions. This is due to the acoustic oscillations which 
become more prominent and their analytic modeling in 
MDM models is more complicated. 



A similar situation occurs also for the dependence of 
TMDAi(k) on h. Since the /i-dependence is included prop- 
erly in kF and k nr , D(k) does not require any correction 
in the asymptotic regimes. Only a tiny correction of D{k) 
in the intermediate range, k (0.01 < k < 1) is necessary 
to minimize the residuals. By numerical experiments we 
find that this can be achieved by multiplying ct\, ...a^ by 
the factors Ci(h) which are approximated by second order 
polynomial on h/O.b with coefficients determined by % 2 
minimization (see Appendix). The fractional residuals of 
D{k) for different h are shown in Fig. 9 (top panel), they 
remain stable in the range 0.3 < h < 0.7. But the frac- 
tional residuals of Tmdm(^) slightly grow (about 2-3%, 
bottom Fig. 9) in the range 0.1 < k < 1 when h grows 
from 0.3 to 0.7. This is caused by the fractional residuals 
of TcDM+b(k) (see middle panel). 

3.3. Dependence on N v . 

The dependence of D(k) on the number of massive neu- 
trino species, N v , is taken into account in our analytic ap- 
proximation by the corresponding dependence of the phys- 
ical parameters k nr and &f (see Eq.(6)). It has the correct 
asymptotic behaviour on small and large scales but rather 
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Fig. 6. The fractional residuals of T(k) given by 
(T(fc) - T C MBfast(k))/T C MBfast(k) for different values of Q„. 



large residuals in the intermediate region 0.01 < k < 10. 
Therefore, the coefficients on (i — 1, ...,4) must be cor- 
rected for N v . To achieve this, we multiply each a, by 
a factor Di(N v ) (~ 1) which we determine by x 2 min- 
imization. These factors depend on N v as second order 
polynomials. They are given in the Appendix. In Fig. 10 
we show the fractional residuals of Tmdm(^) for different 
numbers of massive neutrino species, N„, and several val- 
ues of the parameters fl v , fib, h and z. The performance 
for N v = 2, 3 is approximately the same as for N v = 1. 

4. Performance 

The analytic approximation of D(k) proposed here has 
maximal fractional residuals of less than 5% in the range 
0.01 < k < 1. It is oscillating around the exact numerical 
result (see Fig. 4), which essentially reduces the fractional 
residuals of integral quantities like a(R). Indeed, the mean 
square density perturbation smoothed by a top-hat filter 
of radius R 

1 f°° 

a 2 (R) = — / k 2 P(k)W 2 (kR)dk, 
% n Jo 

where W(x) = 3(sinx — xcosx)/x 3 , P(k) = AkTl 1DM (k) 
(Fig. 11) has fractional residuals which are only about half 



Fig. 7. Fractional residuals of the analytic approximation for 
the MDM transfer function at different redshifts. 

the residuals of the transfer function (Fig. 12). To normal- 
ize the power spectrum to the 4-year COBE data we have 
used the fitting formula by Bunn and White 1997. 

The accuracy of <j(R) obtained by our analytic ap- 
proximation is better than 2% for a wide range of fl u for 
fib = 0.06 and h = 0.5. Increasing fib slightly degrades 
the approximation for N v > 1 , but even for a baryon con- 
tent as high as fib ~ 0.2, the fractional residuals of <r(R) 
do not exceed 5%. Changing h in the range 0.3 — 0.7 and 
N v — 1 — 3 do also not reduce the accuracy of <j{R) beyond 
this limit. 

We now evaluate the quality and fitness of our approx- 
imation in the four dimensional space of parameters. We 
see in Fig. 12 that the largest errors of our approximation 
for (j(R) come from scales of <~ 5 — 10h~ 1 Mpc 4 . Since these 
scales are used for the evaluation of the density perturba- 
tion amplitude on galaxy cluster scale, it is important to 
know how accurately we reproduce them. The quantity 

4 Actually, the oscillations in the error of T(k) are somewhat 
misleading: they are mainly due to baryonic oscillations in the 
numerical T(k) entering the denominator for the error esti- 
mate, so that a slight shift of the phase enhances the error 
artificially. This is why we concentrate on the error of o(R) 
(otherwise the error estimate of T(k) should be averaged, see 
e.g. EH2). 
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as = a(8h 1 Mpc) is actually the most often used value 
to test models. We calculate it for the set of parameters 
0.05 < fl u < 0.5, 0.06 < fl b < 0.3, 0.3 < h < 0.7 and 
N v — 1, 2, 3 by means of our analytic approximation 
and numerically. The relative deviations of erg calculated 
with our TMDM{k) from the numerical results are shown 
in Fig.13-15. 

As one can see from Fig. 13, for 0.3 < h < 0.7 and 
£lh/i 2 < 0.15 the largest error in as for models with one 
sort of massive neutrinos N v = 1 does not exceed 4.5% for 
£l v < 0.5. Thus, for values of h which are followed by direct 
measurements of the Hubble constant, the range of ilb h 2 
where the analytic approximation is very accurate for 
fl u < 0.5 is six times as wide as the range given by nucle- 
osynthesis constraints, {^h 2 < 0.024, Tytler et al. 1996). 
This is important if one wants to determine cosmologi- 
cal parameters by the minimization of the difference be- 
tween the observed and predicted characteristics of the 
large scale structure of the Universe. 

For models with more than one species of massive neu- 
trinos of equal mass (iV„ = 2,3), the accuracy of our an- 
alytic approximation is slightly worse (Fig. 14, 15). But 
even for extremely high values of parameters fib = 0.3, 
h = 0.7, N v = 3 the error in cr 8 does not exceed 6%. 




1E-3 0.01 0.1 1 10 



k (h/Mpc) 

Fig. 9. Fractional residuals of the analytic approximation of 
D(k), TcDM+t(k) and T M DM(k) transfer function for different 
values of the Hubble parameter, h. 

In rcdshift space the accuracy of our analytic approx- 
imation is stable and quite high for rcdshifts of up to 
z = 30. 

5. Comparison with other analytic approximations 

We now compare the accuracy of our analytic approx- 
imation with those cited in the introduction. For com- 
parison with Fig. 12 the fractional residuals of <j(R) cal- 
culated with the analytic approximation of T MDM {k) by 
EH2 are presented in Fig. 16. Their approximation is only 
slightly less accurate (~ 3%) at scales > lOMpc/h. In 
Fig. 17 the fractional residuals of the EH2 approxima- 
tion of TMDAf(k) are shown for the same cosmological 
parameters as in Fig. 16. For fl v = 0.5 (which is not 
shown) the deviation from the numerical result is > 50% 
at k > l/i _1 Mpc, and the EH2 approximation completely 
breaks down in this region of parameter space. 

The analog to Fig. 13 (as) for the fitting formula of 
EH2 is shown in Fig. 18 for different values of fib, f2„ and 
h. Our analytic approximation of TMDAi(k) is more ac- 
curate than EH2 in the range 0.3 < il„ < 0.5 for all Of, 
(< 0.3). For Vl v < 0.3 the accuracies of as are comparable. 

To compare the accuracy of the analytic approxima- 
tions for TMDAi(k) given by Holtzman 1989, Pogosyan & 
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Fig. 10. Fractional residuals for the analytic approximation of 
the MDM transfer function of density perturbations for differ- 
ent numbers of massive neutrino species, N v , in models with 
different values of Q v , fl b and h. 

Starobinsky 1995, Ma 1996, EH2 with the one presented 
here, we determine the transfer functions for the fixed set 
of parameters (Cl v — 0.3, fib = 0.1, N v = 1, h = 0.5) 
for which all of them arc reasonably accurate. Their de- 
viations (in %) from the numerical transfer function are 
shown in Fig. 19. The deviation of the variance of den- 
sity fluctuations for different smoothing scales from the 
numerical result is shown in Fig. 20. Clearly, our analytic 
approximation of Tmdm (k) opens the possibility to de- 
termine the spectrum and its momenta more accurate in 
wider range of scales and parameters. 

6. Conclusions 

We propose an analytic approximation for the linear power 
spectrum of density perturbations in MDM models based 
on a correction of the approximation by EH1 for CDM 
plus baryons. Our formula is more accurate than previous 
ones (Pogosyan & Starobinsky 1995, Ma 1996, EH2) for 
matter dominated Universes (£Im = 1) in a wide range of 
parameters: < Q. v < 0.5, < fl b < 0.3, 0.3 < h < 0.7 
and N v < 3. For models with one, two or three flavors of 
massive neutrinos (N v = 1, 2, 3) it is significantly more 
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Fig. 11. The variance of density perturbations smoothed by a 
top-hat filter of radius R for different MDM models. Numeri- 
cal transfer functions Tmdm are shown as thick lines and the 
analytic approximations Tmdm(^) = TcDM+bD (k) are thin 
lines which are superimposed on the numerical results (The 
residuals are not visible on this scale.). 
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Fig. 12. The fractional residuals of <j(R) defined by 
(a(R)—acMBfa.st(R))/<JCMBfa.st(R) for different values of f2„ 
with all other parameters fixed. 

accurate than the approximation by EH2 and has a rela- 
tive error < 6% in a wider range for Vt v (see Figs. 13, 18). 

The analytic formula given in this paper provides an 
essential tool for testing a broad class of MDM models 
by comparison with different observations like the galaxy 
power spectrum, cluster abundances and evolution, clus- 
tering properties of Ly-a lines etc. Results of such an anal- 
ysis are presented elsewhere. 
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Fig. 13. Deviation of as of our analytic approximation for 
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Our analytic approximation for Tmdm(&) is avail- 
able in the form of a FORTRAN code and can be re- 
quested at novos@astro.franko.lviv.ua or copied from 
http: / / mykonos.unige.ch/^durrer/ 
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Fig. 17. Fractional residuals of the analytic approximation by 
EH2 of the MDM transfer function for the same parameters as 
in Fig.16. For comparison see Fig. 6. 
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Fig. 16. Fractional residuals of a(R) calculated by the analytic 
approximation of EH2 for the same parameters as in Fig. 12. 
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Fig. 19. Fractional residuals of different analytic approxima- 
tions for the MDM transfer function at z = for one flavor of 
massive neutrinos. 



Fig. 18. Deviation of as as obtained by the fitting formula of 
EH2 from numerical results for different values of Q„, Qb and 
h (N v = 1). See Fig. 13 for comparison. 













.„•''' our 




E&H 1997b 




Ma 1996 




P&S 1995 


.... 1 


Holtzman, 1989 

i ■ ■ ■ 



R TH (h'Mpc) 



Fig. 20. Fractional residuals of a(R) calculated with the same 
analytic approximations as Fig. 19. 
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Appendix 

The best fit coefficients aj(f2„), &j(f2„), Cj(f2„), Bi(Q b ), Ci(h) and D i (N l ,): 

ai = 1.24198 - 3.887870 + 28.335920 2 _ 70.9063O 2 + 84.158330* - 41.166670^ , 
a 2 = 0.7295 - 3.61760„ + 21.458340 2 _ 54.63036^ + 70.802740* - 35.209050^ , 
a 3 = 0.28283 - 0. 539870,, + 5.80084O 2 _ 14.182210^ + 16.85506O 4 , - 8.776430^ , 
a 4 = 0.48431 + 1.89092O„ - 4.04224O 2 + 8.096690^ - 10.05315O 4 , + 5.344050^ . 

&i = 0.2667 - 1.670,, + 3.560 2 - 3.10^ , 
d = 0.008 - 0.055O„ + 0.135O 2 _ 0.1240^ ; 

b 2 = 0.226 - 0.47O„ + 1.270 2 _ 1.390^ , 
c 2 = 0.004 - 0.026O„ + 0.053O 2 _ 0.0390^ , 
6 3 = 0.076 - 0.258O„ + 0.2150,2 , 
c 3 = 0.0026 - 0.0205O„ + 0.055O 2 - 0.0510^ , 

b 4 = 0.0158 - 0.055O„ + 0.0228O 2 , 
c 4 = 0.00094 - 0.0072O„ + 0.018O 2 - 0.0160^ . 

Si (fi b ) = 1.202 - 0.2065(O b /0.06) + 0.005(O 6 /0.06) 2 , 
B 2 (Q b ) = 1.033 - 0.036(O b /0.06) + 0.003(O b /0.06) 2 , 
B 3 (Q b ) = 1.166 - 0.17(0/0.06) + 0.005(O/0.06) 2 , 
B 4 (Q b ) = 0.97985 + 0.01525(0/0.06) + 0.00626(O 6 /0.06) 2 . 

d(h) = 1.09 - 0.09(/i/0.5) , 
C 2 (h) = 1.65 - 0.88(/i/0.5) + 0.23(/i/0-5) 2 , 
C 3 (h) = 1.055 - 0.055(/i/0.5) , 
C 4 (h) = 0.981 + 0.03(/i/0-5) - 0.012(/i/0.5) 2 . 

Di(N„) = 1.315 - 0.4317V,, + 0.1167V 2 , 
D 2 {N V ) = 1.108 - 0.2257V„ + 0.1177V 2 , 
£> 3 (7V„) = 1.316 - 0.4327V„ + 0.1167V 2 , 
£> 4 (7V„) = 1.256 - 0.3027V„ + 0.0467V 2 . 
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